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Abstract 

We apply stochastic quantization method to real symmetric matrix models for the 
second quantization of non-orientable loops in both discretized and continuum levels. 
The stochastic process defined by the Langevin equation in loop space describes the 
time evolution of the non-orientable loops defined on non-orientable 2D surfaces. The 
corresponding Fokker-Planck hamiltonian deduces a non-orientable string field theory at 
the continuum limit. 
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String field theory [Q is believed to be the most promising approach to investigate 
non-perturbative effect in string theories. Recently, non-critical string field theories have 
been proposed for c = 0(010 and for < c < 1 |||]. Among these, some 00] 
are based on the transfer-matrix formalism in dynamical triangulation of random 
surfaces 0. While the other P] is derived by using stochastic quantization method |^. 

In the approach by stochastic quantization, introducing the loop variable trM" ^ for 
hermitian matrix models and interpreting the fictitious time as a time coordinate, Jevicki 
and Rodrigues 0] showed that the Fokker-Planck hamiltonian ( or loop space hamilto- 
nian ), in which the loop variable and the conjugate momentum are identified with the 
creation operator and the annihilation operator respectively, realizes the string field the- 
ories which are equivalent to the field theory derived by Ishibashi-Kawai 0. Inspired by 
the work in Ref. 0, we apply stochastic quantization method to real symmetric matrix 



models [|T0| and show that it leads to a field theory of non-orientable (non-critical) strings. 
The stochastic process defined by the Langevin equation in loop space describes the time 
evolution of the non-orientable loops on non-orientable 2D surfaces. The correspond- 
ing Fokker-Planck hamiltonian is a loop space hamiltonian of non-orientable string field 
theories. At the equilibrium limit, it deduces the Virasoro constraint equation for the 
probability distribution functional. The continuum limit of the field theory of discretized 
non-orientable loops is taken for the simplest one-matrix case ( c = ) and deduces the 
continuum field theory of non-orientable strings. 

Let us start with the Langevin equation for one matrix model, 

AM..M = -^^(M),,(r)Ar + Ae.,(r), 
S{M) = -Y, -^A^-"/2trM"+2 ^ 

(1) 

Mij denotes a real symmetric matrix. The fictitious time r is discretized with the unit 
time step At. We consider the discretized version of time evolution Mij{T + At) = 
Mij{T) + AMjj(r) , with the Langevin equation for convenience of stochastic calculus 
and for understanding the corresponding stochastic process precisely. The discretized 



fictitious time development with At precisely corresponds to the one step deformation 
in dynamical triangulation in random surfaces. In the following argument, the specific 
form of the action of the matrix model is not relevant. The correlation of the white noise 
A^jj is defined by 

< A^i^{T)A^ki{r) >5= AT{5u5jk + SikSji) . (2) 

It is uniquely determined^ from the requirement that (1) is transformed covariantly pre- 
serving the white noise correlation (2) invariant under the transformation M — ^ UMU~^ , 
where U denotes orthogonal matrices for the real symmetric matrix models. 

The basic field variables are loop variables (pn = tT{M^)N~^~^ . Following to Ito's 



stochastic calculus [|TT1, we evaluate 



A0n = (pnir + At) - (j)n{T) 



n-2 



-| n—z 

ntr(AMM"-^)A^-^-5 + -^ E tr(AMM'=AMM"-'=-2)Ar-i-f + 0{At^/^) 

^ k=0 



(3) 



The terms in R.H.S. should be of the order At, thus we obtain 
n ""^ 1 

A0„ = Ar-{ Y^ (j)k4>n-k~2 + {n- l)—(f)n-2} + At nY^ QaK+a + AC„-i , 
^ fc=0 ^^ a=Q 

AC„-i = ntr(A^M'^~i)A^-^-t . 



The correlation of the new noise variables appeared in (4) is given by 



(4) 



2 
< ACm-i(r)AC„-i(r) >^= AT—nm < (t)m+n-2{r) >^, (5) 

The new noise is not a simple white noise but includes the value of the loop variable itself. 
In a practical sense, it might be tedious to generate the noise variable. We notice that 
0m+n-2(T) in R.H.S. of eq.(5) does not include the white noise at r but the series of noises 
up to the one step (fictitious time unit At) before. This means that the expectation value 



^ For an hcrmitian matrix Alij in (1), the nose correlation is < A^y (r)A^fc;(r) >{= 2ATSiiSjk ■ 



in R.H.S. should be defined with respect to the white noise correlation up to the fictitious 
time r — Ar. 

We also notice < AC„,(r) >^= by means of Ito's stochastic calculus. In the context 
of SQM approach, the property of the noise yields the Schwinger-Dyson equation by as- 
suming the existence of the equilibrium limit at the infinite fictitious time, or equivalently, 
lim^^oo < A0„(r) >^= . We have, 

< 2 51 4>k4>n-k~2 + ^{n- 1)— 0„_2 +nY^ ga4>n+a >5= . (6) 

fc=0 «=0 

The order of the noise correlation (5), 1/N'^, realizes the factorization condition in the 



large N limit |12|. Therefore we obtain the S-D equation at large N limit for discretized 



non-orientable strings. 

1 n-2 

- ^ < 0fc >5< 0n-fc-2 >i+J29a< 4>n+a >€= . (7) 

^ fc=0 a=0 

This shows that the S-D equation for non-orientable strings takes the same form as that 
for orientable strings at large N limit. The correspondence at the large N limit is exact 
if we define the corresponding hermitian matrix model by replacing all the couplings , 
ga — > 2ga in (1). As a consequence, the disc amplitude in non-orientable strings is exactly 
the same as that in orientable strings. 

The geometrical meaning of the stochastic process described by the Langevin equation 
(4) is the following. The one step fictitious time evolution of a discretized loop, 0„(r) — >• 
0n(''") + A0„(r) , generates the splitting of the loop into two smaller pieces, 0^ and 0n-fc-2- 
The process is described by the first term in R.H.S. of (4). In a field theoretical sense, 
it is interpreted as the annihilation of the loop 0„ and the simultaneous pair creation of 
loops, (pk and 0n-fc-2- The first term in R.H.S. of (4) preserves the orientation of these 
loops, while the second term, which is the characteristic term of the order of -^ for non- 
orientable interaction, does not preserve the orientation. Since the new noise variables 
in (5), A^„_i's, are translated to "annihilation" operators in the corresponding Fokker- 
Planck hamiltonian, the factor 2 in the correlation (5) for the new noise variables comes 
from the sum of the orientation preserving and non-preserving merging interactions. 



Namely, (5) describes the simultaneous annihilation of two loops 0m and 0„ and the 
creation of a loop 0m+n-2- The geometrical picture allows us to identify the power "n" 
of matrices in (pn to the length of the discretized non-orient able loop 0n- We notice that, 
in each time step, the interaction process decreases the discretized loop length by the 
unit "2" . The process which comes from the original action of matrix models extends 
the length of discretized loops. These features are equivalent to the transfer-matrix 
formalism [^] in dynamical triangulation of 2D random surfaces in which the one step 
deformation of a specified loop on a triangulated surface defines a discretized (proper) 
time evolution of the loop. 

The definition of the F-P hamiltonian operator gives the precise definition of a field 
theory of second quantized non-orientable strings. In terms of the expectation value of 
an observable 0(0), a function of 0n's, the F-P hamiltonian operator Hpp is defined by, 

< 0(O)|e--^--O(0)|O >=< 0{ct^^{r)) >g . (8) 

In R.H.S., 0^(t) denotes the solution of the Langevin equation (4) with the initial con- 
figuration 0(0) 7^ 0. The time evolution of R.H.S. is given by, 

< AO(0(r)) >5 = < E 9„^O(0(r)) A0m + ^ ^ a^9„O(0(r)) A0mA0„ >^ +0{^t^'^) , 

m ■^ m,n 

= -Ar < if^p(r)O(0(r)) >5 , 

(9) 

where 9„ = g|-. By substituting the Langevin equation (4) and the noise correlation (5) 
into (9), we obtain 

n>0 

1 1 ""^ 11 

-^n - JJ2YI m(t)rn+n-2T^rri + T^Yl 0fc0n-fc-2 + T^iP' - l)^0n-2 + Yl 9a(pa+2 
^^ m ^ fc=0 ^ ^^ a=0 

(10) 

where 7r„ = -^ . To define the operator formalism corresponding to eq.(8), we introduce 
0m and T^ra as the creation and the annihilation operators for the loop with the length n. 



respectively. Then we assume the commutation relation [vTm, 0n] = ^mn , and the existence 
of the vacuum, |0 >, with 7rm|0 >=< 0\4>m = for m > . In the representation, 
< Q\ =< 0\e^mQ"^*^ and \Q >= UmS{(f)m - <5m)|0 > , the F-P hamiltonian operator 
Hpp in (8) is given by replacing (p^ —>' 4>m, and vr^ —>■ Tim in Hpp in (10) with the same 
operator ordering. 

From the equality (8), the probability distribution function P(0, r) , which is defined 
by < 0{(P{t)) >(.= jYind(j)nO{4>)P{(p,T) , is given by, 

P(0,r)=<0(O)|e--^--|0> . (11) 

The initial distribution, P(</), 0) = Iirn5{(j)m — </>m(0)) , represents the initial value of the 
solution of the Langevin equation (4). Eq.(ll) follows the Fokker-Planck equation for 
the probability distribution, 

AP(0, r) = +Ar^^pP(0, r) , (12) 

where Hpp is the adjoint of Hpp in (10), 

Hpp{t) = -rnin Y^ Xn , 

n>0 

1 1 "~^ 11 

^n= - -J-^'YrmTm4>m+n-2 + 9 XI 4>k4>n-k-2 + ^('^ " '^)j-;4>n-2 

+ Y^ 9a<Pa+2 ■ 

(13) 

In the context of stochastic quantization, the F-P hamiltonian (10) in loop space was 
found for hermitian matrix models. The remarkable observation was that it includes 
the Virasoro constraint [Q. Since the fictitious time evolution is generated by the noise 
essentially, the emergence of Virasoro constraint is traced to the noise correlations in 
eq.(5) which are equivalent to the insertion of matrices into the loop variable, M — »• 
M + ArM™"-^ , in 0„. It generates the transformation [— ArLm-2,0n] = n/S.T(j)m+n-2 , 
which corresponds to the noise correlation (5). In fact, for real symmetric matrix models 
( non-orientable strings ), Ln = —N'^Xn+2 satisfies the Virasoro algebra without central 



extension, 

[Lm, Ln] = (m - n)Lm+n- (14) 

We notice that, although Hpp is not hermitian, one can define an hermitian F-P hamil- 
tonian from Hpp as a direct consequence of the fact that Ito's stochastic calculus auto- 
matically picks up the Jacobian factor [|T^ ( or more precisely, invariant measure ) which 
is induced in the space of loop variables by the change of dynamical variables from a 
matrix to loop variables 0. 

It is also worthwhile to note that the F-P equation (12) realizes the Virasoro constraint 
for the probability distribution. Namely, I/„ = N'^Xn+2 also satisfies the Virasoro algebra 
without central extension (14). Therefore, the F-P equation deduces a constraint equation 
for the distribution function even at the discretized level, justifying the generation of the 
partition function which satisfies the Virasoro constraint at the infinite fictitious time. 



L„ lim P(0, r) = , for n = -1, 0, 1, ... . (15) 

T^OO 

For hermitian matrix models, the Virasoro constraint for the partition function (15) was 



found as the S-D equation ||I^. In the continuum limit, it deduces the continuum version 
of the Virasoro constraints [1^. The expressions (8) and (11) also give a constraint on 
possible initial condition dependence of the expectation value and the partition function 
at the infinite fictitious time limit, such as, limT-^oo -f^Fp[0(O), g^Qy]-P[0, t] = . This 
implies that these quantities may have the initial value dependence up to the solution of 
the Virasoro constraint. 

Now we take the continuum limit. First we introduce a length scale "a " and define 
the physical length of the loop created by 0„ with / = na. Then we may redefine field 
variables and the fictitious time at the continuum limit as follows. 

Gst = N-'a-"" , 
dr = a-2+^/2^r, 
$(/) ^ a-^/20„, 

n(/) = a-i+^/V„, 



(16) 
where we would like to keep the string coupling constant, Gst, finite at the double scaling 



limit 16 . For the existence of the smooth limit from the discretized fictitious time 



evolution to the "continuum" one, we require the condition, ^ — 2 > {). The scaling 
dimensions of all the quantities in (16) have been determined except the scaling dimension 
of the string coupling constant, D, by assuming p|§], 

ArHpp = drTipp , 

m)Mi)] = Ki-i')- 

(17) 
Then we obtain the continuum F-P hamiltonian, 7Y_fp, from Hpp at the continuum hmit, 

^non-or. ^ _ / dl{2Gst / dl'^{l + l')l'U{l')m{l) + / dl'^l - l')^{l')m{l) 

2 JO Jo Jo 

+V^/<i>(/)m(/)+p(/)n(/)}, 

(18) 

for non-orientable strings. By the redefinition (16), the F-P hamiltonian (18) is uniquely 
fixed at the continuum limit except the cosmological term. To specify the explicit form 
of the cosmological term p{l) in (18), we have to carefully evaluate the contribution 
which comes from the 3-point splitting interaction term and the matrix model potential, 
a^^^ I]a=o 9a<Pn+a , at the double scaling limit of the real symmetric matrix model. Here 
we remember that the S-D equation for non-orientable string at large N limit takes 
exactly the same form as the orientable one under the suitable choice of the matrix 
model coupling constants. Since the continuum limit is taken by using the universal part 
of the disc amplitude , we naively expect p{l) takes the same form as that for orientable 
strings. 

To show explicitly this is indeed the case, we take the continuum limit in the real sym- 
metric matrix model with the same procedure in Ref . . We consider the simplest matrix 
potential given by, go = —1/2, gi = g/2,g2 = gs = ... = in (1), which corresponds to 



c = 0. Let us introduce the string field variable 



N z- MN^y^ 
Adz) ^ ^z-i-"AC„. 

n 

To take the continuum limit, we redefine the field variable, 

0(z) = -{z-gz^) + Coz:^a^/^^u) , 
Adz) ^ coz-^a'lHau) , 



(19) 



(20) 

where we have introduced the "renormalized" parameters, z = ZcQ""^ , and g = gcQ~'^^°' * , 
and the "continuum" fictitious time dr = cqz~'^q^I'^ Ar , where Zc = 33(32 + 1) and 
Qc = ^"T are the critical values and t denotes the cosmoloeical constant. The constants 
Co and ci are chosen for convenience. The scaling dimension of all the quantities have 
been determined so that the string coupling, 1/Gst = c^N^a^ , is fixed at the double 
scaling limit [|16|. Then we have the following Langevin equation in continuum limit. 



d^u) = dCiu) - ^drdumuY - -^^du{^[u)) - a-%^iz - gz'f + gz)} , 
= dCiu) - ^drid^^uf - p{u) - ^^dlmu))} + 0{a) , 

<dau)dau')> = ^dA'{^mu)-<i>{u'))}+o{a), 

i\'^a° u — u 

p{u) = 3u^ -— . 

(21) 

We have only picked up the terms in the noise correlation which survive at the continuum 
limit. We notice, since d^{z — gz"^) = 0{a?), the field redefinition in (20) is irrelevant in 
the term of the order l/(A^a^/^). p{u) in (21) is precisely the cosmological term appeared 
in the orientable string. 

By the Laplace transformation. 



roo 

$(m) = / dle~''^^(l) , (22) 

Jo 



/o 
8 



we obtain the Langevin equation which is equivalent to the continuum F-P hamiltonian 

(18), 

d^l) = dcii) + ldT{i f di'^i')<^{i - /') + p{i) + \[G^ti^^{i)] , 

2 Jo * 

< dC{l)dC{l') > = 2dTGJl' < $(/ + /')>, 

p{l) = 36"{l)-^6{l), 

(23) 

for non-orientable string. It is consistent with the naive continuum hmit of its discretized 
version (4) except the term p{l). By the same procedure, the Langevin equation for 
orientable string is given by 

rf$(/) = dC{l) + dT{l f rf/'$(/')$(/ - I') + pil)} , 

Jo 
< dC{l)dC{l') > = 2dTGstll' < $(/ + /')>, 

(24) 

in the hermitian matrix model with replacing the coupling constants, Qa -^ 2ga, in (1). 
As we have shown explicitly, the cosmological term p{l) takes the same form both for 
orientable and non-orientable strings. The field theory of non-orientable strings is also 
consistent with ref. [Q in transfer matrix formalism. The double scaling limit of the real 
symmetric matrix model has been studied in a quartic potential [|I^, while our result 
shows that it happens in the cubic potential as well. We notice that the continuum F-P 
hamiltonian includes the continuum Virasoro generator C{1), 

m = -{ Hdi'^i + i')i'n{i') + -^ tdi'^i - /')$(/') 

Jo 2Gst Jo 

1 ,4(0 +^4^1}. 



2\JGgt 2Gst I 

(25) 

These generators satisfy the continuum Virasoro algebra, [C{1), C{1')] = {I — l')C{l + /') . 
Let us briefly comment on the multi-matrix model cases. We may start from a set of 
Langevin equations. 

AM(p),,(r) = -^^5(M),,(r)Ar + Ae.,(r), 



(26) 

The index p specifies the p-th matrix M{p)ij and the white noise A$,{p)ij. p runs from 1 
to A^O) namely we consider Nq matrices. The size of all the matrices is assumed to be A^. 
The basic loop variable is typically of the form, 

tr(n„M(j9„)"^«) =tr(M(pi)"''iM(p2)"''2...M(j9„)"^"...) , (27) 

The Langevin equation for these variables takes the form similar to the one matrix model 
case. Namely, it includes only linear terms and bilinear terms of loop variables ( anni- 
hilation of a loop and simultaneous creation of a pair of loops ) . The correlation of the 
noise variables is given by a linear combination of loops. It realizes the process where 
various pairs of loops are annihilated simultaneously and a loop is created. In this case, 
the colored noise correlation is still equivalent to the insertion of matrices. Therefore, 
if we identify the length of the loops by the power of the matrices included in the loop 
variable, we may conclude that there exists Virasoro constraint with respect to the loop 
length indices in the loop variables ( string fields ). 

In conclusion, we have shown that the Langevin equation ( or equivalently the cor- 
responding F-P hamiltonian ) for real symmetric matrix models written by the loop 
variables defines the time evolution in the ( non-critical ) non-orientable string field the- 
ory at both discretized and continuum levels. The partition function for non-orientable 
strings satisfies the Virasoro constraint at the equilibrium limit in both discretized and 
continuum level. In the stochastic quantization view point, since the fictitious time scal- 
ing dimension is given by y — 2 = | > for c = 0, we expect that the discretized version 
of the loop space Langevin equation for real symmetric matrix models may provide a 
possible method for numerical calculation of non-orientable 2D random surfaces to sum 
up the topologies of surfaces. 

Note 

10 



After the completion of this work, we found the works, refs. []T7| |T^ [|T^ , where the 



second quantization of master fields are discussed by stochastic quantization method. 
The expectation value of master fields corresponds to eq. (9) in the present scheme. 
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